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I.  INTRODUCTION 

Asynchronous  wavelength  division  multiplexing  (WDM)  systems  have  been 
increasingly  proposed  as  an  attractive  alternative  to  coherent  optical  frequency  division 
multiplexing  (FDM)  systems  [Ref.  1-5].  Although  asynchronous  WDM  systems  with 
direct  detection  do  not  have  the  channel  capacity  of  coherent  optical  FDM  systems,  they 
are  much  less  costly  to  implement.  Furthermore,  present  filter  technology  enables  the 
designers  to  tightly  pack  the  channel,  resulting  in  asynchronous  dense  WDM  systems  that 
can  provide  aggregate  bit  rates  of  many  terabits  per  second  (1  Tb/s  =  1012  b/s). 
Asynchronous  dense  WDM  systems  are  particularly  attractive  in  the  area  of  undersea 
surveillance  where  hundreds  of  sensors  and  data  collection  sites  are  envisioned  being 
merged  onto  single-fiber  superhighways  through  massive  data  fusion.  Other  applications 
call  for  relatively  low-bandwidth  data  collection  over  many  months  to  be  dumped  quickly 
to  a  remote  recording  site  in  a  matter  of  minutes.  This  "collection-and-dumped" 
compression  can  demand  total  data  rates  on  the  order  of  hundreds  of  Gb/s.  Long  distance 
between  data  collection  sites  and  a  remote  recording  site  requires  the  use  of  optical 
amplifiers.  Therefore,  it  is  necessary  to  pack  all  channels  within  the  optical  amplifier 
bandwidth. 

An  asynchronous  dense  WDM  receiver  with  on-off-keying  (OOK)  modulation  can 
be  modeled  as  shown  in  Fig.  1.  Conceptually,  the  analysis  involves  two  main  operations: 


1)  a  convolution  operation  to  evaluate  the  signal  at  the  output  of  the  optical  filter,  a 
Fabry-Perot  (FP)  filter  in  our  investigation,  and  2)  the  integration  of  the  output  of  the 
photodetector.  Evaluation  of  bit  error  probability  by  the  numerical  analysis  of  these  two 
operations  has  been  carried  out  in  [Ref  6],  with  a  number  of  approximations  made  to 
reduce  the  computational  complexity.  In  this  investigation  the  FP  filter  is  shown  to  be 
well  approximated  by  an  RC  filter  within  the  frequency  range  \f-f0  \  <  FSR I  20tt,  where 
FSR  is  the  free  spectral  range  of  the  FP  filter  [Ref.  7]  and  /0  is  the  FP  filter  center 
frequency.  For  example,  given  FSR  =  3800  GHz,  the  approximation  works  very  well  for 
|  f-f0\<  60.5  GHz;  that  is,  the  effects  of  adjacent  channels  within  a  121  GHz  bandwidth 
centered  at/0  must  be  included,  while  all  others  can  be  neglected.  This  simple  model 
agrees  well  with  [Ref.  6]  as  demonstrated  in  Section  HI.  Furthermore,  this  model  enables 
us  to  obtain  a  closed  form  analytical  expression  for  the  bit  error  probability  for  which 
numerical  results  can  be  obtained  with  little  effort.  Our  investigation  shows  that  this 
simple  model  provides  accurate  results  as  compared  to  those  in  [Ref.  6]  for  bit  rates  up  to 
2.5  Gb/s  when  the  effects  of  two  adjacent  channels  are  included  with  FSR  =  3800  GHz. 
Actually,  when  the  power  penalty  relative  to  single  channel  operation  is  1  dB  or  less,  there 
is  virtually  no  difference  in  the  effect  of  four  or  two  adjacent  channels.  Thus,  for  this 
power  penalty  criterion,  this  simple  model  can  handle  bit  rates  up  to  5  Gb/s  for  a  FP  filter's 
FSR  =  3800  GHz. 

In  Section  II  the  closed  form  expression  for  the  decision  variable,  and  consequently, 
the  bit  error  probability  assuming  all  channels  are  bit  asynchronous  as  in  [Ref.  6]  is 


derived.  Section  III  presents  the  numerical  results  which  include  the  bit  error  probability 
versus  the  signal-to-noise  ratio  as  function  of  the  FP  filter  bandwidth  and  channel  spacing, 
and  the  power  penalty  (relative  to  single  channel  operation  without  filtering  or  with 
filtering  but  no  intersymbol  interference)  versus  the  channel  spacing  as  a  function  of  the 
bandwidth.  Finally,  a  summary  of  results  appears  in  Section  rV. 


H.  ANALYSIS 

The  receiver  model  for  the  asynchronous  dense  WDM  system  is  shown  in  Fig.  1. 
The  desired  signal  is  filtered  by  a  Fabry-Perot  (FP)  filter  that  rejects  adjacent  channels. 
The  photodetector  is  assumed  to  have  a  responsivity  (H  (AAV).  The  detected  current  is 
amplified  by  a  low  noise  amplifier  that  contributes  a  postdetection  thermal  noise  n{i)  with 
spectral  density  N0  (A2/Hz).  The  decision  variable  at  the  output  of  the  integration  is 
compared  to  a  threshold  a  to  determine  whether  a  bit  zero  or  bit  one  was  present. 
A.      INPUT  SIGNAL 

For  convenience,  we  designate  channel  0  as  the  desired  channel,  and  channel  k  as 
an  adjacent  channel  where  k  =  -Mil,  ...,  -1,  1,  ...,  M/l  with  M  an  even  integer.  We 
consider  the  equivalent  lowpass  (complex  envelope)  data  signal  in  channel  0  and  channel  k 
as  follows: 

b0(t)  =     i    bojPiit-il)  (i) 

bk{t)  =    I   bkje^-^Prit-ilT+Xk))  (2) 

l=-L 

where 

T:  bit  duration 

b0ie  {0,1 }:  bit  in  channel  0  in  the  time  interval  (/T,  (/+1)7) 


bkl  e  {0,  e*k)  is  the  Ith  bit  in  channel  k  in  the  time  interval  (IT+  xk,  (/  +  1)7  +  xk) 
(J)t :  a  phase  offset  between  channel  k  and  channel  0  and  is  assumed  to  be  uniformly 

distributed  in  (0,  2n)  radians 

©k :  radian  frequency  spacing  between  channel  k  and  channel  0  with  cok  =  -co.k 
Xk  :  is  the  time  delay  between  channel  k  and  channel  0  and  is  assumed  to  be 
uniformly  distributed  in  (0,T). 

The  function  PT(t  -  iT)  is  defined  as 


FT{t-ll)-{0       otherwise  (3) 


In  both  (1)  and  (2),  the  non-negative  integers  L0  and  L  represent  the  number  of  bits  in 
channel  0  and  k,  respectively,  that  proceed  the  detected  bits  b0Q.  The  received 
asynchronous  dense  WDM  equivalent  lowpass  signal  at  the  input  of  the  FP  filter  is  given 

by 

Mil        ,— 

r(t)=JPb0(t)+    S     JPbk(t)  (4) 

k=-M/2 
where  P  is  the  received  optical  power. 


B.       FABRY-PEROT  FILTERED  OUTPUT  SIGNAL 

The  FP  filter  can  be  characterized  by  the  following  equivalent  lowpass  transfer 
function  [Ref.  1,7] 


pe-j2T<f/FSR  l_p 

1-p  1-^-p 


• 


l-pcos(||)+ypsin(||)         l  P 


(5) 


where  p  is  the  power  reflectivity,  A  is  the  power  absorption  loss  (zero  for  an  ideal  FP 
filter)  and  FSR  is  the  free  spectral  range.  For  \f\<  FSR/20iz  and  assuming  A  =  0,  we  can 
approximate  H  (/)  as  follow: 


('-PH/lf      l+Jl$k 


\/\<FSR/20n  (6a) 


where 


C  = r-^  (6b) 


The  free  spectral  range  FSR  can  be  related  to  the  full  width  at  half  maximum  (FWHM) 
bandwidth  B  and  the  finesse  F  of  the  FP  filter  as 


nJpB 

FSR  =  -r—  =  BF  (?) 

i-p  v  ' 


Thus  if  the  signal  is  bandlimited  to  |/|  <  FSRJ20n,  we  can  truly  approximate  (5)  with  a 
single-pole  RC  filter  with  the  following  transfer  function  and  impulse  response 


i+rr 


h(t)  =  ce~ct,   t>0  (9) 

Figures  2a-b  show  the  magnitude  and  phase  (radians)  of  H{f)  of  the  FP  filter  in  (5)  and 
its  single-pole  RC  filter  approximation  given  in  (8)  for  p  =  0.99,  F  =  312.6,  B  =  12.16 
GHz  and  FSR  =  3800  GHz.  Note  that  as  the  frequency  increases,  the  phases  of  the  FP 
filter  and  the  RC  filter  differ  markedly,  but  the  magnitudes  of  their  transfer  functions 
remain  identical  and  attenuate  rapidly.  When  |  /  |  >  FSR/20K,  the  magnitude  of //(/)  is 
very  small,  and  therefore,  the  effect  of  adjacent  channel  interference  beyond  this  frequency 
range  is  negligible.  Fig.  3  shows  the  normalized  impulse  response  of  both  FP  and 
single-pole  RC  filters.  In  summary,  the  above  approximation  is  valid  for  asynchronous 
dense  WDM  analysis  when  the  filter  finess  F  is  large  or  equivalently  the  FWHM 
bandwidth  B  is  small  since  the  equivalent  lowpass  signal  must  be  bandlimited  to  about 
|  /  |  <  FSR/20n. 

This  approximation  has  been  used  in  [Ref  5]  to  study  spectral  efficiency  of  optical 
FDM/ASK  systems,  which  involves  the  evaluation  of  the  decision  variable  for  worst-case 


analysis  using  the  eye  diagram  technique.   Since  we  are  interested  in  the  detected  bit  b  00 
in    the  time  interval  (  0,T),  we  consider  the  output  filtered  signal  s(t),  0  <  t  <  T  given  by 

s(t)  =  sB(.f)+sISi(t)+sAci(f),    0<t<T  (io) 

where 

sB  (/):  desired  signal 

s]SJ(t ):  intersymbol  interference  (ISI)  signal 

sACI  (t):  adjacent  channel  interference  (ACI)  signal 
These  signals  are  evaluated  using  (4)  and  (9)  as  follows  (detailed  derivations  in  Appendix 
A): 

o 
=  JPb0,o(l-e-ct),     0<t<T  (11) 


_!  (H-D7* 

sisi(t)  =  JP    2    b0 j     J     h(t-X)dk 

i=-Lo  iT 


=  JPe-0'    X    b0J(e(m)cT-eicT),  0<t<T   (12) 

i=-L0 


Saci(0 

_     MR.         -2  <;+I)7-«t 

=  J?    I,    [Z,  btj      J      /z(?-?.)e^(X-Tt)^  0<a<r 


t=-M/2   i=-l  IT+Tk 


t 


+  64,-1    J  Kt-xye^'^dk  0<t<tk 

-T+xk 

+  Vi    ?    Kt-K)e>**^dk  Tk<t<T 

~T+Tk 

+  bk/\h(t-\)e^-^dk]  Xk<t<T 


JF  *f   [£  ti(0,  7)6*  /  e"c'  ^  e{o&*W«*)(p<rt*kyr  _  ^ 

k=-M/2    l=-L  '  1+7-r 

/fc*0 


+7i(0,  x*)6*-i  e"c'  ^  (e^*)'  -  e^*)*-^*)) 


+71(1*,  7)6*_i  e'cr  £Z!^  e(<*fl>k)^(\  -  eH°*»ky!) 


+it(Tt,  7)A*,0  e"c'  ^  (e(c+>*)<  -  e(<*»*)**)]         03) 


0<t<T 


where    71  (/1/2)  —  (o  otherwise     's  a  u™1  Pu'se  between  /,    and  t2 


The  FP  filtered  output  s(t)  is  detected  by  the  photodetector  which  produces  a 
current  of  d  \  s(t)  |2  Amps.     This  current  plus  additive  white  postdetection  thermal    noise 
current  from  the  amplifier  is  integrated  by  the  integrator  to  obtain  a  decision  variable  for 
the  threshold  detector. 
C.      DECISION  VARIABLES 

The  decision  variable  Y  appearing  at  the  integrator  output  consists  of  the  signal 
component  X  and  noise  component  N 


where 


Y  =  X+N  (14) 


T 

X=\     As(f)Vdt  (15) 

0 


T 

N=\n(t)dt  (16) 

o 

We  note  that  N  is  a  zero  mean  Gaussian  random  variable  with  variance  N0T.   Substituting 

(10)-(13)  into  (15)  we  obtain  the  signal  component  X  as  a  function  of  the  three  parameters 
Tfc,  cT,  and  o\r,  which  represent  the  effect  of  intersymbol  interference  and  adjacent 
channel  interference. 

X=      £  PTb20o  [1  -jf{\  -e^r)  +  ^(l  -e~lcT)\ 
+  ££L(\-e-2<T)[  t   boAe(M)cT-ewT)]2 


10 


+  ^(1  -e~2^  I     I      I  Jw^L^»^-  l)l2 

icy  fc=-M2  /=-/_       '      !+/-# 


/fc*0 


oT  k=-M12   m=-MI2         (l+Z-rXW^)         Lg/   *  '*  co^co^ 


pT      M/2        M/2      bk-ib^.y  e-XWk-""**")     («,-(<**»*)**_  i)  e(c*jakX-T^k) 

+  £^  k=M,2    m=M/2         (l^)(l-7^)  [  H^¥ 

(e-(c-Ju>k}ts_i\e(c-ja>myrT+tm)  (e-2c*5_l)  g  («*/«»  jfeX"^*)  e(c-ytom)(-7,+xw) 

1   i  c 


pT        M/2  M/2        bk-Xb*     .   (e-2«g-e-2^  g-Xto^-tomx^) 

+  £EI?        2  Z  "^    ST ^ 

cT  k=-MI2  m=-M/2                   2(l+y-^)(l-7T) 

fc*0  m*0 


r^(c+7©fc)T,fc/j  _g-(c+7ojk)7,\e(c-7a>l„)Tlw/j  _  g-(c-7G>«)7Y| 


pr        M2  M/2        6*n6*  A  e~j(<»k*k-<*mzm)         _ 

k=-MI2    m=-M/2        (1+J— )(!-./—)  ,    --  w„ ©**©« 


pT      M/2        M/2       bk^b*m0e-Jia^k-^^m)        ^j0kyzk  (e-ic+j<»m)T_  e-i<*i*my*g) 
*  k=-M/2    m=-M/2         (H/¥xw¥)  H/¥ 

+  wf  l  ] 


11 


a?2PTd      f      ¥2  %2      bk^b^Qe-^k^m){l_e^«0T)e{c^k),k 

+  "?  —zrRe  {     L,  L,      —h — 

cT         k=-MH  m=-M/2                   (i+hrXW5?) 

btO  rmtO 


e-(c+J<t>m)xg_  e-(c+J(om)T  (e~2c?g _e~2cTj  e(c-j«>m)xm 


->DT,               M2  M2       -2    birib*    ,  e'^W*"1^  (e(c+J<*^T-\)  e(c+jmkX1T*k) 

+  <f^Re{    I        Z      2  -^^ J     . 

cT            k=-MI2  m=-MJ2l=-L                            fl+hrXl-/?) 

k*0  m*0 


[      f^  2  ]} 


+  ^  ^S/te  {    S        IS  - ^r — ^ 

cr  fc=-M2    m=-M/2l=-L  (l+Jir)(HT) 

k*Q  m*0 


bicjb*m_l(e-2"m-e-2cT)  (l-e"<^->'")r)  e(Ha>mytm 

I  9 


bk,ib*m0(e-2czm-e-2cT)  e^J'om^m 

?  ~  J  / 
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+  £  —fRe  {     Lt  L 


cT'  k=-MI2    m=-M/2  (1+7^X1-;^) 

fc*0  m*Q 

bic-\b*    ,(e~<-c~jak^m-e~('c~jai^kVl-e~^c~J<iiw)T)  e(^c~J<i>m^m 


[■ 


h4 


bk-\b*    Ae~2cXm-e~2cx^  (\-e-<c-J(»»>)T)  e(c+Jak^~T^0  e{c-ju>m)T-m 


bK-lb*mQ{e}{(*k-<»>n*k-e>(,»k-<»>n)-"n) 

J         c 


bk-ib*m  Q(e~{p~iai<y*n  -e^&kfttyic-jomiytm 


1-/* 


bk-\b*  Je<^mm^m-e'^c+jam)'lk)e^c+J(i'k^-'rj^k) 


i+;— 


(The  above  whole  term  will  be  zero  if  ©£  =  CO/w). 


+  *  ^o,o(i-^    ^    b,Ml+m  -  ewT) 


cT 
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+  d—^rKe  {     la         Li 


cT  k=-M/2   l=-L  2(1+7-/) 

toO 

k=-M/2         l+j-r  J-r  H~r  2 

toO 


ifc=-M2  l+/-r  2 

toO 

NP-      bofibw  e~JmFk  ^eJOkT-eJo^k       e-(c-J(i>k^k-e-<^-JVik^T 
£=-M2         l+j-r  J-r  H~r 

toO 


(2e_ct't-2e-c7'-e-2c^+e-2cr)e(c+-'a,A:)x't 

~~ 7  J    J 


+  ^^    S    6o,-(e('+1)cr-e/cr) 

/=-Z-o 

n     ,     M2  =2     ^(l-*-™)    (g(^M)g(^)/rg^ 

ite  {     2-        2    — ST 

k=-M/2  l=-L                            2(l+j-f) 

toO 

^  Bl  0)7,  ^ 

£=-M/2  l+j-r        l-j-r  l 

toO 


] 


^W  ^Ar.O      rg-cxA:-e--/(0^g-(c--'<0/t)7'  e-"k-e"k  e-2cT-i    X 

toO 


where  t5  =  min  (t*,  xm)  xg  =  max(T*,  im) 
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D.       BIT  ERROR  PROBABILITY 

For  a  detection  threshold  a  and  an  ISI/ACI  bit  pattern  b  =  {&u  ,  b0i  };  /  =  -L0,  ...-1; 
/  =  -Z,  ...,  0;  k=  -M/2,  ...,  M/2,  fe  0;  the  conditional  bit  error  probability  of  the  OOK 
signal  represented  by  the  Gaussian  random  variable  Fin  (14)-(16)  is  given  by  [Ref.  8] 


where  Q  (X)  is  defined  as 


Q(x)  =  -p=r  J  e~y2/2dy 

^y   '       flZ    x  *  (19) 


and  X0  and  A',  are  the  values  of  X  in  (17)  for  b00  =  0  and  b00  =  1,  respectively.  The 
average  bit  error  probability  P,  is  obtained  by  taking  the  expected  value  of  Pe  (b)  in  (18) 
over  all  bit  patterns  b.  The  minimum  bit  error  probability  is  obtained  by  optimizing  over 
the  threshold  a.  In  summary,  given  a,  we  calculate  the  following  expectations: 


Pe  =E  {Pe(b)}  (20) 

b 


Pe,mm=mmE{Pe(b)} 
a    b 


=  mjn  -^tr  _  2      p(b) 

a      l  u    2M patterns 
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P(b) =k-f-  )Mci)w 


2V  In  J     yT 


In        T 

[  J  •  •  J    J  •  •  •  J  Q( 7= )d§-MI2  ■  ■  •  «q>M/2   dX-M/2 . ..dlMIl 

MO     MO  v^o/ 


MO     A/  0  V  W 


(21) 
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m.  NUMERICAL  RESULTS 


In  this  section  we  present  the  numerical  results  for  a)    bit  error  probability  versus 


signal-to-noise  ratio  Z  =  RPJT/No  as  a  function  of  the  normalized  channel  spacing 
(normalized  to  the  bit  rate) 


7  =  g  =  (A/)7-      (22) 


where  Af=  G\/  2nk  is  the  equal  channel  spacing  in  Hz,  and  b)    power  penalty    versus 
normalized  channel  spacing  /as  a  function  of  the  FP  filter    parameter  cT=(n/Jp)BT 
where  B  is  the  filter  full  width  at  half  maximum  bandwidth(FWHM).    For  our  model  to  be 

valid  we  set   M  =  2,    that  is,    we  constraint  the  signal  to  be  bandlimited  to  within  two 

adjacent  channels,  thus  we  incorporate  the  degradation    Hv  the  two  nearest  adjacent 

channels.    We  observe  that  for  bit  rates  of  1  Gb/s  or  less,  oui  model  is  valid  up  to  M=  10 

and  very  little  difference  is  observed  between  M=  2  and  M—  10.  Also,  we  observe  that 
there  is  little  difference  between  M=2  and  M-  4  when  I>10  for  bit  rates  up  to  3  Gb/s. 
In  all  results  we  set  L0  =  2  and  L  =  1 . 

Figures  4-5  show  the  minimum  bit  error  probability  versus  signal-to-noise  ratio 

(  Z  )  for  cT  =  5  and  10,  respectively.  Along  with  each  curve,  we  also  show  that  of  the 
synchronous  case,  i,e,  Zk  =  0,L0  =  2  and  L  =  0,  and  that  of  a  single  channel  (SC) 
operation  without  filtering  or  with  filtering  but  without  ISI.   In  Fig.  4  we  observe  that   a 

large  degradation  occurs  due  to  ISI  for  cT=  5  which  represents  a  narrowband  filter.    As 
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the  FP  filter  bandwidth  is  made  larger  as  in  Fig.  5  with  cT=  10,  the  IS1  is  reduced  but  the 

ACI  increases. 

In  our  model,  we  are  constrained  to  M  =  2  for  the  case  under  consideration.    We 

use  an  FP  filter  with  FWHM  5=12.16  GHz,  free  spectral  range  FSR  =  3800  GHz,  finess 
F=7typ7(l  -p)  =  FSR/B  =  312.6,  and  c  =  38.4  GHz;  1/7  =  2.56  Gb/s,  and  cT  =  15.  If 
the  power  penalty  related  to  single  channel  operation  is  to  be  less  than  1  dB,  then  by 

Fig  6,   a  minimum  normalized  channel  spacing  of  I  =  12  must  be  used.     Equivalently 

from  [Eq.  (22)],  we  could  use   Af  =  I/T  =  12*2.56  GHz  =  30.4  GHz,  (i.e.  the  channel 

spacing  is  a  multiple  of  the  bit  rate).  In  this  model  the  farthest  adajcent  channel  for  M=  2 

is  twice  the  channel  spacing  which  is  60.8  GHz.   This  verifies  the  assumption  I  f-f0  I  < 

FSR/20iz  =  60.5  GHz,  where  f0  is  the  FP  filter  center  frequency.  This  result  agrees  well 

with  that  in  [Ref.  6;  Figs.  6,9,  M/F  =  0.4,  a  =  0,2].  Thus  we  incorporate  the  degradation 

caused  by  the  two  nearest  adjacent  channels.  We  observe  that  for  bit  rates  of  1  Gb/s  or 

less,  our  model  is  valid  up  to  M  =  10,  and  very  little  difference  is  observed  between  M  - 

2  and  M  =  10.  Also,  we  observe  that  there  is  little  difference  between  M  =  2  and  M  =  4 
when  I  >  10  for  bit  rates  up  to  3  Gb/s.  In  all  results  we  set  L0  =  2  and  L  =  I. 

Figure  6  also  shows  the  power  penalty  for  an  asynchronous  dense  WDM  system 

relative  to  a  single  channel  operation  at  the  minimum  bit  error  probability  of  10"15.  This  is 

the  required  additional  signal  power  (dBW)  for  the  asynchronous  dense  WDM  system  to 

be  able  to  operate  at  the  10'15  bit  error  probability  achieved  in  the  single  channel  system 

with  a  SNR=l2dB.  The  asynchronous  dense  WDM  system  is  ZS7-limited  at  2. 15  dB,  0.95 

dB,  0.6  dB,  and  0.5  dB  in  power  penalty  for  cT  =  5,  10,  15,  and  20,  respectively.    It  is 
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seen  that  for  a  2.3  dB  power  penalty,  the  normalized  channel  spacing  can  be  as  close  as  / 

=  6  (i.e.,  a  channel  spacing  of  six  times  the  bit  rate)  for  cT  =  5.  If  the  power  penalty 
criterion  is  1  dB,  the  normalized  channel  spacing  is/«  12  for  cT=  10,  15,  20.  We  remark 
that  although  the  exact  transfer  functions  of  the  FP  filter  is  used  in  [Ref.  6],  a  number  of 

approximation  have  been  made  to  obtain  numerical  results.  The  approximations  are  1) 
the  ISI  is  obtained  by  modeling  FP  filter  as  a  single-pole  RC  filter  [Ref.  6,  Eqs.  (4)  and 
(36)],  2)  approximating  the  finite  integration  with  an  infinite  integration  in  the  calculation 
ofACI  [Ref.  6,  Eq.  (15)],  and  3)  the  beat  interference  is  ignored.  On  the  other  hand,  the 
ISI  and  ACI  in  our  investigation  are  obtained  by  modeling  the  FP  filter  as  a  single-pole  RC 
filter,  using  finite  integration  and  including  the  beat  interference.  Since  the  results  in  our 
investigation  and  in  [Ref.  6]  agree  well,  we  conclude  that  approximations  are  quite  valid. 
We  also  note  that  our  results  also  agree  well  with  the  simulation  carried  out  in  [Ref.  1, 
Fig.  17]. 

The  above  numerical  results  shown  in  Figs.  4-6  are  obtained  with  an  optimized 
threshold  setting.  Figure  7  shows  the  power  penalty  for  fixed  threshold  a  =  £PTI2  which 
is  the  same  optimum  threshold  for  single  channel  operation  (midpoint  between  the 
received  power  for  bit  zero  and  bit  one).  It  is  seen  that  the  performance  of  an 
asynchronous  dense  WDM  system  is  quite  sensitive  to  a  for  a  narrow  band  filter.  An 
additional  1.5  dB  is  observed  for  cT  =  5  for  I  >  8,  and  0.4  dB  for  cT  =  10  for  /  >  12. 
Negligible  degradation  is  observed  for  cT=  15,  20  for  />  16. 
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Figures  8-9  show  the  power  penalty  versus  normalized  channel  spacing  as  a  function 

of  FP  filter  parameter  cT  for  the  worst-case  analysis  with  optimal  threshold  and  fixed 

threshold,  respectively.     The  worst-case  bit  pattern  is  fixed  to  produce  the  minimum    X] 

and  maximum  X0  where  Xx  andX0  are  the  values  of  X   in  Appendix- A  equation  (4-21) 

with  b0  0  =  1  and  b00  =  0  respectively. 

We  observe  that  the  power  penalty  for  the  worst-case  analysis  is  only  slightly  larger 

than  that  of  the  exact  analysis  for  /  >  10  shown  in  Fig.  6.   Similarly  the  power  penalty  for 

the  worst-case  analysis  with  fixed  threshold  is  only  slightly  larger  than  that  of  the  exact 

analysis  with  fixed  threshold  for  /  >  10  shown  is  Fig.  7.    The  reason  for  this  is  that  for 

large  channel  spacing  (/  >  10),  the  ACI  effect  is  small,  so  the  ACI  bit  pattern  has  a  small 

influence  on  the  power  penalty. 

Figure  10  shows  the  normalized  optimal  threshold  for  the  exact  analysis  shown  in 
Fig.  6.  It  is  observed  that  a  «0.4  for  /  >  10.  Note  that  the  normalized  optimal  threshold 
for  the  single  channel  operation  is  a  =  0.5. 
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rv.  CONCLUSIONS 

We  have  presented  a  simple  model  for  the  analysis  of  asynchronous  dense  WDM 

systems  employing  an  external  OOK  modulator.    The  only  approximation  that  we  use 

involves  the  modeling  of  the  Fabry-Perot  filter  by  a  single-pole  RC  filter  assuming  the 
equivalent  lowpass  signal  is  bandlimited  to  the  frequency  range  l/l  <FSR/20n.  This 
model  enables  us  to  obtain  a  closed  form  expression  for  the  bit  error  probability  which 

previously  could  only  be  obtained  via  numerical  analysis  [Ref.  6].  For  an  FP  filter  with  an 

FSR  around  3800   GHz,  our  model  can  include  the  ACI  effects  of  two  adjacent  channels 

for  bit  rates  up  to  2.5  Gb/s.   Our  numerical  results  show  that  this  model  agrees  well  with 

that  in  [Ref.  6]. 
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Figure  1  :OOK  receiver  structure 
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Figure  2:  Spectral  characteristics  of  the  Fabry-Perot  filter  and  the  approximated 
single-pole  RC  lowpass  filter 
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Figure  3:  Normalized  rimpulse  response  of  the  Fabry-Perot  filter  and  the  approximated 
single-pole  RC  lowpass  filter 


24 


10* 


ASYNC 


SYNC 


16 


Figure  4:  Probability  of  bit  error  versus  signal-to-noise  ratio  as  a  function  of  normalized 
channel  spacing  for  cT=5 
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Figure  5:  Probability  of  bit  error  versus  signal-to-noise  ratio  as  a  function  of  normalized 
channel  spacing  for  cT=\0 
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Figure  6:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  cT 
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Figure  7:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  cT  with  a  fixed  threshold  a=0.5 
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Figure  8:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
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Figure  9:  Power  penalty  versus  normalized  channel  spacing  as  a  function  of  Fabry-Perot 
filter  parameter  c^for  worst-case  with  fixed  threshold  a=0.5 
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Figure  10:  The  normalized  optimal  threshold  versus  normalized  channel  spacing  as  a 
function  of  Fabry-Perot  filter  parameter  cT 
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APPENDIX  A 


Derivation  of  Formula 


/.      INTRODUCTION 

Desired  channel  Channel  0. 

Adjacent  channel   :  Channel  k. 

k  =  -M/2,  ...  -1,  1,  ...  M/2,      M:  even. 

Bit  in  channel  0  in  ith  time  interval  (/T,  (i+l)T)  :  b0i  e  {0,  1 } 

Detected  bit  in  channel  0  in  (0,T) :  b00  6  {0, 1 } 

Bit  in  channel  k  in  lth  time  interval    (IT+  Tk,  (/  +  1)7>  xk)  :  bk,ieJ(*k(t~Zk) 

bkJe  {0,e>*k},j=J=T 

($k     frequency  spacing  between  channel  k  and  channel  0,  CO>t  —  —  ®-k. 

fok  ;    phase  offset  between  channel  k  and  channel  0,  uniformly 
distributed  between  (0,  In). 

Tic',     time  delay  between  channel  k  and  channel  0,  0<t^<T. 


o 
Data  signal  in  channel  0  :  bo(t)  -  £    ^o,/  Prif  -  iT)  (1-1) 

i=-Lo 
0 

Data  signal  in  channel  k :  bk(t)  =  £  ft*  iei^(t-xk)  p^  _  qT+  ^  (1_2) 

L0,  L  :  positive  integers, 
p  rA  =  |     i,o</<r 

*  7w        l      Otherwise 

The  received  signal  at  the  input  of  the  FP  filter  of  channel  0  is  : 
r(t)=JPb0(t)+    E     JPbk(t) 

k=-MI2 

P  :  received  optical  power. 
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The  output  of  the  FP  filter  is: 


oo 


r0(f)=  J  h(t-X)r(X)dX 


—00 


oc  °o 

Mil 


i —    r  / —     Mil       r 

=  Jp  J  h(t-X)bo(X)dX+JP     £      J  h(t-X)bk(X)dX 


-oo  k=—M/2   -oo 

too 


00 


=  VP6o,oJ  h(t-X)PT(X)dX 


—00 

-1  °° 


+  V?   Z   6o.i  I  h(t-X)PT(X-iT)dX 

i=-Lo  -oo 

^_     M2  -2  °f  ,,       , 

+VF    S       {  I  &*,/   J  h(t-X)e^x^PT(X-(lT+Tk))dX 


k=-MI2  l=-L  -oo 

too 


00 


+&*,_!  J  h{t-X)e^k^-^PT{X-{-T+xk))dX 

—00 

oo 

+6*,o  j  h(t  -  X)e^x-^PT(X  -  xk)dX  }  (i 


-oo 


h(t)  is  the  equivalent  lowpass  impulse  response  of  the  FP  filter  of  channel  0. 
Since  the  detection  interval  is  0<t<T,  we  only  need  to  evaluate 

Sft)  =  r0  ft),  0<t<T  and  write    Sft)  as  : 

Sft)  =  SB(t)  +  Slsl(t)  +  SAC1  ft)  o<t<T 

SB  ft)  :  Desired  signal  in  channel  0,  the  1st  term  (i=0)  of  (1-3) 
Sisi(t)  :  Inter-symbol  interference,  the  2nd  term  ft  <  -I)  of  (1-3) 
SACIft):  Adjacent  channel  interference,  the  3rd  term  (Z)  of  (1-3) 

k 

All  evaluated  in  0<t<T. 
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//.     CALCULATE  S&  (t)  and SM  (t) 

In  this  section  we  are  to  calculate  the  easy  part  ofSft)  —  SB  ft)  and  SISI  ft). 

A.      FP  filter : 

The  lowpass  equivalent  transfer  function  of  FP  filter  is 

1-p  \-A-p  1-p  l-A-p 


H(0  = 


1  -^-J2n/IFSR    '       1_D  '  .2:1/.     .       .    .271/.    *       1_D 

1   ?e  P  l-pcos(^)+/psin(T^)  P 


FSR'  J^        ^FSR- 


P 

A 

FSR 


power  reflectivity. 

power  absorption  loss  (  0  for  ideal  filter). 

free  spectral  range. 


Since  f«  fsr  (for  operating  frequency  range),  we  can  approximate 
H (f)  as  (assume  A  =  0) : 


1-p  1-p  i 

VS  271/"  271/p  271/" 

(H>>t/Pi5      l+J7Tp^R-      l+J- 

FSR(l-p) 
where  c  =  — p — 

for  FP  filter,  we  also  have  : 

FSR       "^  Jp 

-g-  =  -^-  ,  B  :  full  width  at  half  maximum  bandwidth  (FWHM) 
of  half  power  bandwidth. 

Using  inverse  Fourier  transform,  we  find 

ce~c\  t  >0 


hft)={ 


0,  otherwise 
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B.      Calculate  SB  (t) : 

oo 

SB(t)  =JPb0,o  J  h(t-X)PT(X)dX 


—00 

t 


=  JPbo,o\h{t-X)dX 


o 

t 


=  JPb0,o\ce-c«-X)dX 


o 


=  JP  b o,o  ce~ct  \ecXdX 


o 


=  y?6o.o(i-o  0</<r 


(a)  For  Bit  "1"    :     SB]  (t)  =  JP  (I  -  e~ct)      0<t<T 

(b)  For  Bit  "0"    :    SB0(t)  =  0  0<t<T 


C.      Calculate  SISI(t) 

SIS1  (t)=JP    2   b0,j  J  h(t-X)  PT  (X  -  iT)  dX 


-isi- 

-l  °° 


i=-L 


o  -00 


_1  (!+l)T 

=  JP    Z   *o,i    J    Kt-X)dX 

i=-Lo  iT 

-i         (HF 

=  JP    I,  b0,i    J    ce-^-^dk 

i=-L0  iT 

-1 

=  JPe~ct   E   b0ti(e^+l^T-eicT)dk  0<f<r 
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///.    CALCULATE  S^ft) 

In  this  section  we  are  to  find  SAC1  with  the  procedure 

A.  When  /  <  -2 

B.  When  /  =  -  1 

C.  When  /  =  0 

D.  Summary  for  SACI  (t) 

A.      When  I  <  -2 


Integrate  the  first  term  of  bk  (t)  in  (1-3),  we  have  : 

Mil      -2  °° 

JP     Z     Z  bkJ  \  h(t  -  X)  e^x-^PT(X  -  (JT+  xk))dX 


fctO 


_2  (l+l)T+Tk 

=  JP  Z  Z  **,/     J     *(*  -  ^)  e^^dk 
k    l      '    zr+c* 
_2      (/+i)r+** 
=  V?  Z  Z  6*,/     J      ce"c^  e*°>&-x*>dk 

k     l  IT+u 

_2  (l+l)T+xk 

=  V?  Z  Z  hj     J     ce^**"***  e<rt*&dk 
k    l      '    /r+r4 

=  Vp  Z  Z  **•»  e^^*T^e^*x^*)(e(^*)r_ 1) 

k     I    l+J-r 

=  #ZE^  e-('-*)(e(<+>*)(/^  _  e(c4/G,,)/r)  0<r<r 
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B.      when  /  =  -1 


B-l     /  =  -L  0<t<T 


Integrate  the  second  term  of  bk  (t)  in  (1-3)  : 

Mil  °° 


yfP     ti    btr-i  j  h(t  -  X)  e^^Prfr  -  (-T+  ik))dk 


k=-M/2 

t 


=  JP2,bt-i    J    h(t-X)e^x-x^dX 


k  -T+z 

t 


=  JP  Z  h-i   J    ce<ct+^kXk)  e(c+^k)XdX 


-T+xl 


£  1+/  c  J 


£-2     /  =  -l,  xk<t<T 

Mil  °2 


k=-M/2 

Ik 


=  JPYibk,-\    J    h(t-X)e^x-^dX  I 

=  JplL^k  e-«-^  (l-  *-<•*»*>*)  t*  <  *  <  r    1 
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C.      when  1  =  0 

C-l.  I  =  0,  0<t<T> 


JP  X  bk,o  j  h(t  -  X)e^k(k-Xk)  PT(k-Tk)  =  0  0<t<  xk         (3-4) 


C-2.  l  =  0,xk<t<T 

Integrate  the  third  term  of  bk  ft)  in  (1-3) : 

oo 

JPYibkfi  \  h(t-X)eJ^x-Xk)  PT(X-xk)cfk 

"•  -oo 

t 

=  V^  X  6*o  f  h{t-X)e^k(K-x^  dk 
*      '  I 

k  ^  Xk<t<T  (3-5) 

=  Jp^  -5*£-  (eju>k(t-iic)  _e-c{t-xk)^  Xk<t<T  (3-5a) 

k    1+^ 


D.      Summary  for  SAa  ft) 

Define  :  K(t2 ,  t2)  =  Utl  ft)  -  Ut2  ft),  f  >  0,  h<t2 

where     Uc(t)={\2  ^  >  0 

is  the  unit  step  function. 
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men   71  (tp  t2)  -    \$  otherwise 


is  a  unit  pulse  between  tl  and  t2 . 

hence  for  0  <t  <T        SAC1  (t)    can  be  written  as 


—      Mil       -2 

^c7  rC  =      JP     X      X  71(0,  7)ijy  e-°l  ^  eW<okW«k){e(c+j<ok)T  _  j j  (3 

fcrt 

(3-6) 

M/2  ,ffl 

+  Vp     X     n(0,  T*)6*,-i  e"*  ^£  (g(^ycD,)/  _  gN^K^))        (3. 

/fc=-M/2  1+^~ 

M/2 

+  [  V?     X     nfo,  7)6*,-i  ^  ^  eW&Hl  -  e-+*^T)         (3- 

yfc=-M/2  l+J~ 

fctO 

M/2  T 

+  V?     X     niXk^b^e^^ie^^-e^^^)]  (3-i 

>fc=-A//2  1+>~ 

fcrt) 


—  VT  -L-  7Xk  _i_  T^  /_  tT 


-  Z0  +  Z0  +  Z%k{-  ZTjk_,  +  ZXk0) 


4-> 

(3-1) 

2?-* 

(3-2) 

<->  -» 

(3-3) 

4,  -» 

(3-5) 
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IV.     WORST-CASE  ANALYSIS 


A.  Worst-case  for  S(t). 

B.  Worst-case  for  £  jQ  \  S(t)  \2  dt  -  output  of  integrator. 

C.  Special  worst-case  analysis. 

D.  Special  worst-case  at  the  output  of  the  integrator. 


A.      Worst-case  for  S(t) 

A-l     Worst-case,  ISI :    boi  =  b_,     £<,  =  °° 

From  (2-2),  we  have  : 

S^,{t)  =  JPb-e-c,{ecT-\)  S  e*T 

h  =  -i 


ecT  I  <  1 


=  Jpb. 

.e^'ie07- 

=  Jpb. 

.e-c'ie07  - 

-l)  e" 

=  Jpb- 

e-*'(ecT  - 

■lK<L 

=  JPb. 

.e-<* 

0<t<T  (4.1} 
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A-2     Worst-case  ACI  :       bkl  =  bk  ;  =  bk0  =  b,     L-^° 


53^0  =  JPb    Z    {7i(0, 7)^  (e^^  - 1)[  X  (e^<°*>r)1 

+V?"  6  X  «0,  ik)-^  (ejm{t-Zk)  -  e-c{t~Xk)  e-{c+^)T)  } 
+JP  b  X  Mxk,  7)—^-  (e-*C-**>  -  e^-^  e^+J^7) 

+  n(Tk,  7)-^  (^'"^  -e^*>)  } 


JPbZ  7C(0, 7)-^-  e-*^  e^^7 
k  1+7-r 


+V?  6  X  71(0,  T*)^sr  (eA»*<^*)  -  e^~^  e-{c^(ak)T) 
+JPb  X  n(xk9  T)-±-  (eJ^*-^  -  e~<^  e-{c+^k)T) 


k  l+J-r 

k  1+i7 


y[Pb^n(0,r) 


eW  0<t<T 
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1 

B.      Worst-case  for  d*  j  IS(r)l2  dt  -  output  ofintesrator 


Output  of  photodetector  :    £  \S(t)\2 

\S(t)\2=  \SB(t)  +  SJt)  +  SACI  (t)  |2 

=  \SB(t)\2  +  \SJt)\2  +  \SACI  (t)  \2  +  2Re{SB  (t)  S'ISI(t)  } 

+  2Re{SrB(t)  SACI(t)  }  +  2Re  { S'ISI(t)  SACI(t)  }  (4-3) 

B-l     From  (2-1) 

l^c  (t)  \2  =  \SB(t)\2  =  Pb„,o  0-2e"  +  e2V,  0<t<T,  (4-4) 

T 

f  \SB(t)\2dt  =  PTb20fi[l  - 1(1  - e-<T)  +  ^(1  - e-^r)]  (4-5) 

0 


B-2     From  (4-1) 


\S^(t)\2  =  Pb2  e2ct,  0<t<T, 

T 


\\S^{t)\2dt=P-£(\-e-™), 


0 
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B-3     From  (4-2) 


\S7C02  =  P\b\ 


eia>lc('-zk^  e-ja>m(t-xm) 


Mil  Mil 

k=-M/2  m=-M/2     ( ]  "tf "e")  ( *  ~^"c") 


k     m 


(i+/?)d-y^) 


0<t<T, 


k    m  (l+7-r)(W^) 


=  { 


ss 


,-^0)^-fflffltffl) 


/t     m  Jy  c 


j&-r)(i+&)(i-m 


(eJ(<dk-®m)T  _   |\         CO*  96  (0TO 


P7lal2S 


*       1+^ 


co*  =  com 
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B-4 


s~-ct       „-2ct\  n  ^  (4-11) 


2Re  {SB(i)  SJS7  (0}  =  2Pb0,ob-(e-ct  -  e'^)  0<t<T 


\2Re  {SB<t)S£®}  dt  =  P-^(l  -e~cTY  (4-12) 


B-5 

2Re  {SB(t)  S7CI(t)}  =  2Pb0,oRe  {b  (1  -  *-**)£  T^  } 

k      l+J— 

=  2Pb0,0Re  {b  E  ^  (e*>*  -  *-<<->*)')}  0  <  t  <  T    (4'13) 

k   l+J— 

A  jt  J  c  J  c  J  c 


B-6 


2Re  {S£®  S7C!(t)}  =  2Pb.Re  {b  Z  ^  *"<**■*}    0  <  /  <  r       M  i« 

A:    i+/—  (4"15) 

2Re  {57,7(0  «)><fc    =  ^  **  <*  2  ^  (1  -  €■*>+#)}       (4.i6) 
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C.      Special  worst-case  analysis 

(£>k  =  2nkI/T,  I  :  positive  integer 

x*  =  0 

b0j  =  b-e  {0,1} 

bkJ  =  bk,o  =  be  {09e*} 

Re{b}e  {0,  cos<|>  } 

\b\2  e  {0,1  } 

£7    J|5b(0I2*  =  ^AJ,0  [1-^(1  -e-cT)  +  ^-T(\-e-2cT)] 


Q2    \\S^{t)\]cdt  =  P-^{\-e-2cT) 


C-3     in  (4-9) 
T 


0 

i 


\\SwAca(i)\ldt  =  PT\b\2Y, 


©7.    -> 


o 


C-4     T 

—    \2Re  {SB(t)S^\t)  }scdt=^-  (1  -e^f 

o  (4-12) 

QS    )lRe  {S*B(t)S7CI(t)  }scdt=2PTb°»;e{b)  (e-«-  1)£  -*-  (4-19) 


_         ±|  2Re  {Sl7(t)S7a(t)  }scdt  =  22^*1(1  _  e-cT)Yj  -^ 
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D.      Special  worst-case  at  the  output  of  the  integrator 

T 

X?CC=^  j\Swc(t)\ldt 

0 

=  £  [  all  terms  in  section  C]. 

+H(l-e-2<r) 

+  *i^(l_e-cr)2 

+  [l*l2+^^^W(*_-*0,o)]2:-^ri}        (4-21) 
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V.      EXACT  ANALYSIS  -GENERAL  CASE 


Now,  we  are  to  encounter  the  most  complicated  situation,  for  clarity,  we 
first  list  the  construction  talbe,  for  which  we  are  to  follow  : 


S(t) 


Sac/0  SJt)  SB(t) 


Z£        Zl"         ZT0  Table  5-1 


z^°        4,-i 


\S(t)\2 

I 

DIC  B 


2Re{S-s,(t)SAa(.l)}     2Re{Sl(l)Sja(l)}        2R*{Sb(i)S;s,M}        \SAa(tP        \S,s,0)\2       1&0)I: 


C-6 C-5 CM C-3 C-2 C- 1 

2Re{Zl>Z?t}  2Re{ZlZVk)    2Re{ZT0zh     IZll2     IzJM2    |z„rl2 


C-3-3 C-3-2 C-3-1 

Table  5-^  2Re{ZiMZl'M}     \zTJ>      \z\J* 

We  first  calculate  each  term  in  Table  (5-2),  then  integrate  it.      For  each  term 
belongs  to  \Saci(0\2,  we  also  discuss  two  cases:  ol>*  =  com  and  u)it  ^  o)OT   . 
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Since  the  computation  is  very  complicated,  here  we  just  write  down 
the  results  for  each  term,  and  omit  the  intermediate  steps. 


A.      J  \SB{t)\2dt  =  PTbloll  - ^p  +  £±p]  (4-5) 


B.      From  (2-2) 


-l 
\SISi{t)\ 2  =  Pe~2ct[  X  bo,i(eii+l)cT ~ eicT)]2  0<t<T  (5-1) 

i=-Lo 


-1 

I 

/=-Z,0 


]\SJS!(t)\2dt  =  %(\-e-2<T)[±  boAem)cT-eicT)]2  (5"2) 


C      Calcualte  \SAC/t)f  -  refer  to  (3-6) 


\Z\V  =Pe-2ctn(0, 7)1  E2**/  ^T  el^W^XeW*&-  l)l2        (5-3) 


r  -2 

\\zl\2dt  =  ^{\-    *T)\X2bkj^eM*W«'XeW^T-l)\2     (5-4) 


48 


C-2  is  =  min  (x*.xm) 


IzJ'l2 


PEE 7C(0, Zs)e-2abK-1b^_l ;^~Z.(e<>*«*  - e«™) 


(i+y-#)(i-7^) 


re(c-jiOm)(  _  e(c-jam)(-T+TmU 


,-j((Ok1k-iiimtm) 


PEE n(0, x^Ai,., ^;,-ie,7;'"";ri[e;(mt-'0"" - e-<c+-"°""  e(^>(-^> 


£     w 


d+y-rKW^) 


_^-(c-yco)k)^(c-7tom)(-7,-Kw)  +e-2cte(e*j(ok)(-T+ik)e(c-jG>m)(-T+zm)-i 

{p-5f 


T 


\\Zx0k\2dt=$\zlk\2dt: 


(a)     k  *  m,  CO/:  *  co 


w 


P7  y  y    ,  ,  *  e-J((Qkik-amXm) 


Hr^T  (e*0*-""*'  - 1) 


"A?-*) 


j !_  (g-fOt/Um)!,  _  J\  e(C+/(Ot)(-r+tt) 


1+/? 


!_  (e-(<;->(Ot)ts  _  j\  e(c-ja„Yr-T-nm) 

'■~J~c~ 

1  (e-2rcs 


|  \  e(c+j(ok)(-T+ik)e(c-j(am)(-T-vcm) j 
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(b)     k-  m,  (Ok  -  o 


m 


_  PT  y    l6*,-i  I2    r  _J_  /.-(c+/co,)t,  _  i  \e(c+j<ok)(-T+ik) 

-  cT  f  1+(?)2  LCT*  +  1+y?  V  l>e 


_I^-2crfc  _  j\  g(c+ya> jkX-r+Tjk)  e(c-;(o/t)(-r+Tfc) 


=  §Sr^[«*-kl-eto*)«-2eJ 


cr  r  i^)2 


+2tfe{— ^  (<r<<;+>>*>7'  -  e^^X-^))}]  (5-6b) 


7-5    3  terms  for  \Z[k\ 

C-3-1     Tg  =  max(T*,  xOT) 


IZ^  I2  =  P  Z  S  7c(xg,  7)  e'^h^b^ 


*       e-A»kxirnm'tm) 


T  m  ""  *  (i+/?)(i-;^) 


[^ (c+> ^)Tfc n  _  g -(c+/co *)7\  e (c-ja m)xmn_e -(c-j(o m)Tyi  /5_7j 


\\Zl_fdt=\\Z\K_?dt 


(a)  k^m,  (Ok  *  co 


m 


PT  y  y  bk.-\b*m-\  ,  _2cTg      ~-2<?7\    g~-/((0^^wT,w) 

crrt    2    l  ;(i+7?)(i-^)' 

[g (c+ja k)*k (\-e -{<»& k)T\  e (c-ja w)t„  / 1  _  e -(c-;co m)7\i  (5-8a) 
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(b)  k  =  m,  (£>k  =  co 


I**, 


m 


C-3-2 


(5-8b) 


\zTJ2=P?lTln(Tg,T)bk,0b*mfi 


k     m 


^-./(co^-co/nTm) 

(1+/— )(l-y— ) 


rej((ok-(om)t  _  e-(c+j(x>m)t\  e(c+j<»k)?k  _  e-4c-jakyt  e(c-j<om)xf 
jr£  -2cte(c+j(i>  k)xk  e  (c-yco  w)t  „  j 


(5-9) 


\\Zl^dt^\\Zl^dt 


(a)  k*m,  (Ok  *  co 


w 


=  77  Z  2  &£,0& 


*  g-Xwjfct/fc-aw*"!) 


,0^^7,0    ,.       WA:, 


•torn- 


[ 


teJ{<»k-<*m)T  _  ej(<*k-<*m)lg\ 


■  £ZZZ  (e-(c+jam)T  _  e-(c+jam)Tg\ 
/^mm)T"  (a-iH^)T  _  p-{c-j<*k)x, 

01 

1-/- 


(« 


) 


_!(£-2cr  _  e-2crg\  g(c+y<Dt)Ttg(c-7<om)T„-| 


(6)  £  =  AW,C0£  =  C0 


/w 


+2^g{e(c^(e"(c+>(0^)y'  -  g-C^"^* 


)}] 
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C-3-3 

2Re{Zl_,ZVK0 }  =  2PRe{I,  I  n(x g.T)b t.,blfif^f=- 


re-ic+j<am)t  g(c+/o*)Tjfc/|  _  e-{c+j(ak)T\ 

-e ~lct  e (c"-/co m)Xm  e (ct/co k)Xk (\-e "(c+;co ;)7)1 }  (5-H) 


\2Re{Z[k_xZl, }dt  =     J  ^{Zi^, }dt       fQT      k^mmdk  =  m 


o  T« 

=  ^*e{I  S  ^.^^r:,  .(1  -«***)  .«*** 


_i^-2cxg  _  ^-2*7)  e(c-j<om)Tm]  j  (5_12) 
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-2 

2Re{ZlZzj  }=2P/te{SS  2  *(0,Tm)MV-i: f^'TT' 

£     w  /=-£  (1+/— )(l-y— ) 

["^ -(c+/<o m)te(c+j<a k)(lT+xk) <e (c+Jco fc)7  _  j  \ 

-£~2c'  ^(c_>'»K-Z'+T'")^(c+/C0't)(/r+T^(^(c+-/C)^)r-  J)]}  /5I 


]2Re{Zlzlk}dt=  J  2Re{Z0Z0k}dt  for  k*m  and  k=m 


=^{ee  2  m;-i  ret;  (*Wm*)r-  d^^ 

fl  _e-(c+/o)/w)xm  (\_p-2cT.m)        .  .    „         . 

rii-f ii-f p(c-ja>m)(-T+zm)-\). 

L        1+y^  2  " 
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C-5 


2Re{ZT0zrk }  =  IPRe^  1 1  n(xm,T)f^=-  (e^^_  l)e^tw^ 


[bk,ib*m-i  e-2c'e(c-Ja''")^(l-e-<c--'a-")T) 

+bklb*m0  e-**** 

-bub'm0  e-toe*?*-*-]}  (5-15) 


\ 2Re{ZTQZrk  }dt  =    j  2Re{ZT0Zrk  }dt 


h    h* 
r     *,/  w.-i   /e-2cxm  _e-2cT\(y  _  e-(c-j(om)T\e(c-j(om)T, 


'■TJ    c 


bJ^±  (e-2cTm  _  e-2cTy{c-janyzm]  j  (5.16) 
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C-6 


2Re{Z0  Z{k  }    The  product  is  only  survived  when  xm  <  t  <  Ty 


e-j((okzk-amzm) 


2Pife{E  E^  (*/»>**) 


-bk-\b*       e-2cte(c+J(°k)(-T+':k)e(c-jGim)xm/Y  _  e~(c-j(om)T\ 

-bK-\ b*m,0  e^c+J*m)te(c+j(0k)(--T+tk) 
+6*,_i Z^ 0  e-2cte(c+JakX-T*k)e(c-JG,m*m]} 


(5-17) 
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\2Re{ZlkZT;k}dt=  \  2Re{ZlkZrk}dt 


2PT 
cT 


/te{2  2 


e/((Oj.iru«t«) 


.bi-\b„ 


r  k-~]   ™-~]  f£-(c-m>)Tm  _  e-{c-j®k)ik\n  _e-(c-Ja>m)T\e(c-j(»m)Tm 

h        h* 

°t,-iyi/ -2crw  -e~lcik)(\  -e~(c~Jam^T)e(c+JCOk^~T+'Zk>}e(c~J(*m>)i:' 

h      h* 
kH_mfi_  ,  j(pk-<nm)xk  _     j(ayk-(£>m)xm\ 

.Ac         c    / 


lit 

i-y-r 


1      ■« 

1+7- 

h       b* 
^  K-i  mfir-2cxm  _e-2cTk\e(c+jak)(-T+Tk)e(c-jam)Tm'i  j 


(5-18) 


D.      J2Re{SB(t)  S*ISI(t)}  dt  =  \lSB(t)  Sisi{f)  dt 


-l 

I 

i=-Lo 


=  2Pb0fi  S  b0,i  (e('+1)cr  -  eicT)  j  (<TC'  -  e"2cr)  rff 


-l 


%  ba,o{\  -  e~cT)2   I  Z,o,,(e('+1^-e'cO 

/=-Z,o 


(5-19) 
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E.      2Re{S*B(t)SAa(t)} 

=  2PRe{m    7i(0,  T)(e~ct  -e-2ct)boobki^^(e^c+J^T  -  \)e^^1T+x^ 

k     I  1+/— 

+Z  7i(0, xk)b0 obk-x'-^ie^  - e-^^1  -  (e~ct  - e-^tyc+J(,k){-T^ 
k  '     l+j-r 

+1 71(1*,  T)b0 obk-ie-^(e-ct  - e~2ct)(\  - e^^e^"^ 

k  '     i+/— 

+Z  n(xk,  T)b0iobk,o  e-^[eJ(i>kt  -  e<°-**#  -  (e~ct  -  e-2ct)e^^}  }       (5-20) 

k  1+7— 


\2Re{S*B(f)  SAa(f)}dt 


1+7- 


CV  ,  ,  2.  l  i  y — £, 


£  1+7—  7—  1-7— 


e->jfe'jfe  r/2e-"*-2e-cr-e"2^+e-2cr 


+i  bofih,-!*-^  v2*  k-^'r  *•* >  (i-  ^<)>(^] 


1+7^ 


^    ,  ,  e"Wfc       /eJa>kT-eJ'»kxk.  ,e-<c-J'»k)xk-e-(r-J0>k)T. 

+2,  t>0,0Ok,0—^  L( — 7^ — )  -  ( — ~k ) 

k  l+J—  J—  1-7— 
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F.      2Re{S]SI(i)SAa(t)} 


=  2P[  Z    boAe(l+l)cT-ewT)] 


j=-L 


o 


-2 


&{ZZ  71(0,  T)bki  e~lct  ^(e^*>r  -  l)e^'^lT^ 


kl=-L  1+/-T 


-J'akxk 


*  i+/- 

+S  7i(Ti,  7)  6*_,  ^  e"2c'(l  -  e^<">*)Vc+><)Ti 

+E  <**,  7)  6*o  ^  («"<»**  -  e-2"e<e+^^)  } 

*  IW  (5-22) 


1 

\2Re{S'lsI(t)SACI(t)}dt 


4rI    6o,,  (e(,+1)cr  -  eicT)Re{YL  h^^p.^j^)T  _  1)e(^tvrg«t 

c      i=-Lo  k  I       2(1+/—) 

ik  i+y—         l-y— 
k  20+/-r) 


+  S  -^sr  Hsr  [e'^  -  e'**''Xke'<c-**$r)  -  \{e~CTk  -  eCXke~2cT)  ]  } 
k   l+J—    H— 


(5-23) 


58 


VI.    EXACT  ANALYSIS  -SPECIAL  CASE 

Now,  we  put  two  constraints  into  all  terms  of  exact  analysis  to  form  a 
special  case  : 

1 .  CD  £  =  ~      I  :  positive  integer. 

2.  Tyfe  =  I    for  all  adjacent  channel    k. 

For  the  purpose  of  normalization  for  future  use,  we  set    a  =  j  ,  or  x^  =  t  =  a  7 


so  that    (DfrT*  =  GDfrT  =  =-jraT  =  2ltkla 


Still,  we  follow  Table   5-2. 


r 


A.      jl5B(Ol2^  =  P^o[l-^p  +  i^P]  (6-1) 


B. 


I5a(/)l2*=#(l  ~^2cr)  [  i  bo.,(e{M)cT-e'cT)  ]2 

7=-Z.o 


C.      Special  case  for  ACI 


C-l 


Q  K  I—  L  J  cj- 

=  <EI{e2«cT  _  e2(a-,)cr)(ecr  _  1)2  |  j|  2gL|  2  (6-3) 
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C-2 


j\r0\2dt=l\zi\2dt 


]m 


bk-\  b*   , 


=  (c7)2pri;S- 

t     m     (cT+j2nkI)(cT-j2nmI) 
l_e-J2ida(k-m) 

*■  j2id{k-m) 

o  ~cT_e  (a-l)cT+j2nmIa 


+ 


+ 


cT+j2nmI 

p  -cT_p  (a- 1  )cT-j2nkIa 

cT-j2nkJ 

e2(a-l)cT_e-2cT 


2cT 


]  (6-4a) 


(by&=  W,CO£  =  CO 


w 


*— i  I/),     J2  /„2(a-l)cr    „-2c7\ 


«  ~         0-cT    0{a-\)cT+flnkIa 
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C-3     3  terms  for    \Z^\2 


C-3-l        J  \z[_x  Vdt  =  J  \z\.x\2dt 


(a)      k^m,  co^co 


w 


*>*.-!  *C.-I 


=  IcTWtl  -  ^(1  -  ^^i:  2  ^^      (6-5a) 


(b)      k  =  m,(dk  =  (D 


m 


=  \cTPT{\  -  e~cT)2(\  -  e2^-^cT)Z 


K-,  I2 


J-  (cT)2+(2nkI)2  (6-5b) 


C-3-2 


J  7 

\\z{Q\2dt  =  \\ZTx()\2dt 


o 
(a)      k^m,  co/:  ^  co 


w 


bk,ob*m0 


-  (Cl)     rl  Zj  Li  (cT+j2iikI)(cT-j2TimI) 


fe-j2Tda(k-m)_]\ 


+ 


+ 


j2nl(k-m) 

(„  (a- 1  )cT+j2nmIa  _  i  \ 

cT+jlnml 

(e  (a-l)cT-j2nkIa  _  j  \ 

cT-j2nkI 


+     2cr      1  <6"6a> 
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(b)         k=m,   (£>k  =  ®m 

=  (cT)2PT  E      '2    '     2  [  1  -  a  +  ( 


^    (c7)2+(27t£/)2   L  2c7 

„(a-l)cr+/2»t/a     i 

+2^  <'   cT^-1}  )  (Mb) 


C-3-3 


J2^e  {z,7;.,^  }df=j2fe  {z7;,^  }A 


(a)      k^m,  g>£  *  co 


w 


2{cT)2PTRe  {£  2  -^%^- 

km   (cT+j2nkI)(cT-A™I) 


(b)     £=/725   C0£  =  CO 


w 


=  2(cT)2PTRe{ZbV%°^ZT) 

y      J  l^     (cT)2+(2nkl)2 


(l_e(a-\)cT+j2iikIa}  g2(a-l)cr_j 

I-        cT+j2%kI  h        2cT        J  *  (6-7b) 
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C-4 


\lRe  {  Zl2$  }dt  =  J2Re  {  ZjZj"  }dt 

o  o 

(a)     k&  m,  (Sik*-  ©m 


2(cT)^T Re  ^zZi^^Tfl 

k    ml    (cT+J2nkI)(cT-J2nmr) 


eacT+j2nmIa_i  e~cT _e{2a-\)cT  . 


(b)      &=/?2,   G)£  =  CDm 


C-5 

T 


\lRe  {  ZT0zr  }dt    J2Re  {  Z\Z^  }dt 


cT+j2nml  2cT         J/  (6-8a) 


eacT+J2nkIa_i  e  -cT  _e(2a-l)cT 

l  cr+/27i^    +        2^        J  '  (6-8b) 


0  T 

v      J  jf  j  {cT+j2nkI){cT-j2nmI)  L  2cT 

+ '  cT+J2nmI +  2^  J>  (6-9) 

C-6 


\2Re{Zz0zr}dt=    \2Re{Zt0Zr}dt=0 


(6-10) 


63 


D.      ]2Re{SB(f)S*lsI(t)  }dt 


o 

=  ^(\-e-cT)2.    Z   boAem)cT-e'cT)  (6-11) 

i=-Lo 


i 

E.     J2Re{SB(t)SACi(t)}dt 


o 


1  *T  —  2(cT+j2nkT) 


+E 


b0,obk-\cT      (\-e-J2nkIa)        (e^T-e^1*1")        ^_e-acT^2e(a-\)cT 

cT+j2nkI    *~~~]2nid        '       cT-jlttkl  2cT         ~ 

b0,obk,-i  (2-2e^-^cT-e-a£T+e^a-2>cT)  (\-e~cT) 


2(cT+j2itkT) 

fC 


+^  cT+j2itkI  *■     ;27c/t/        h        cT-j2nkl  2^T 

(6-12) 

T 

F-    \2Re{S]sl(t)SACi{i)}dt 

o 

=  2/T[  S  Z)0,,(e('+1)cr-e'cr)] 

/=-Z-o 
fV  v  bkJ(l-e-^)  {e^)cT_e<KT)elcT 

{2^y  2(cT+j2nkI) 

+ Y  cT+j2nkI  I-      cT-j2nkl  2cT  -I 


+E 


^(e^-^-^Hl-e"^) 


2(cr+/27t£/) 

y       bkficT         (e-«cT_e-cT-j7xkIa)        ^-acT_e(a-2)cT^  (6"13) 

+ Y  cT+j2nkl  ■■        cT-j2%kI  2cT  J  J 
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APPENDIX  B 


Programs 


%  exact  analysis  with  optimal  threshold  signal  out  of  the  integrator,  the  signal  X  contains  desired 
signal  and  ACI  &  ISI  and  postdetect  noise  ,The  formula  for  the  BER  is  : 


1 


2^^     lUpattern 


2      P(b) 


where  p(b)  =  Hi)MU)M 


2  v  In  '     v  T 


L  J  •  •  •  J    J  •  •  •  J  Q{ J== )"<P -Mil  •  ■  •  d§Mll   ClX -Mil  ■  ■  .dlM/2 

MO     MO  V^o/ 

.     f         f     f         f  nt^l  Q-M2 ■  ■  ■  $M/2  T-M/2 ■  ■  -tjW/2 )-<* x  jX  i,  j  ^7  1 

+  J  —J   J   •  J  &r( Trf^ )«<P-m/2  •  •  .acpM/2  "T-M/2  •  •  -dlMIl  J 

MO     MO  V^°y 

M  is  the  number  of  adjacent  channels.   Here  we  assume  all  channels  have  the  same  phase  (j)  and 

time  delay  I  i.e,  §k  —  §,  t/c  —  T    hence  for  our  model  the  power  of  (1/27C)  and  (1/T)  are  fixed 
to  1,  therefore  we  only  have  two  integrals  and  two  arguments. 

M=4;k=[-M/2:-l  l:M/2]; 
m=k; 

%  produce  the  controlled  matrix  b  to  control  64  different  bit  patterns 
ml=[  zeros(l,32)  ones(l,32)  ]; 
m2=[  zeros(l,16)  ones(l,16)  zeros(l,16)  ones(l,16)]; 
m3=[]; 

fori=l:4 

m3=[m3  [zeros(l,8)  ones(l,8)]]; 

end 
m4=[]; 

fori=l:8 

m4=[m4  [zeros(l,4)  ones(l,4)]]; 

end 
m5=[]; 

fori=l:16 

m5=[m5  [zeros(l,2)  ones(l,2)]]; 

end 
m6=[]; 

fori=l:32 

m6=[m6  zeros(l,l)  ones(l,l)]; 

end 
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b=[ml;m2;m3;m4;m5;m6]'; 

%  signal  to  noise  ratio  range  in  dB 
RPTN_DB=[  10:0.5:25];  %RPTN_dB 
ppp=10.A(0. 1  *RPTN_DB); 
lenl=length(ppp); 

%  solalph  function  is  provided  by  Professor  Randy  L.  Borchardt 

n=10; 

[bpx,wfx]=grule(n);  %bpx=bpy  ,wfx=wfy 

%  single  channel 
BER0=0.5*erfc(ppp/8A0.5); 

pp=[];  thresh=[];     %thresh  is  not  normalized  threshold 

for  CT=[5  10  15  20]     %  cT  is  Fabry-Perot  filter  parameter 
ifCT=5 

qqq=[2:12  17  20]; 
elseifCT=10 

qqq=[3:12  17  20]; 
elseifCT=15 

qqq=[4:12  17  20]; 
elseifCT=20 

qqq=[5:12  17  20]; 
end 


pe=[];  threshl=[]; 

for  I=qqq 

BER=[];  thresh2=[]; 

for  RPTN=ppp 

ap=linspace(0, 1 82, 1 1 ); 
x3=0; 

%  BER  minimun  don't  care  about  the  scale 
fori=l:64    x3=x3+... 

solalph(,xx00',0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+.. 
solalph(,xxll\0,2*pi,2,bpx,wfic,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 
lef=ap(ind)-16.6; 


%  approximated  thresholds 
%  first  few  loops  to  find  out  the  threshold  which  make  the 
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if  lef<0  ,lef=0;  end       %  to  aviod  the  threshold  go  beyond  the  negative  side 

ap=linspace(lef,ap(ind)+ 16.6,1 1); 
%  two  more  time  to  find  alpha 
x3=0; 

for  i=  1:64     x3=x3+... 

solalph(,xx00',0,2*pi,2,bpx,wfx,0, 1 ,2,bpx,  wfx,ap,b,CT,i,I,k,m,RPTN)+. . . 

solalph('xxll\0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 
lef=ap(ind)-3 . 1 ; 

if  lef<0  ,lef=0;  end 

ap=linspace(lef,ap(ind)+3 .1,11); 
%  three  more  time  to  find  alpha 
x3=0; 

for  i=l  :64 
x3=x3+... 

solalph('xx00\0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 

solalph('xxl  l',0,2*pi,2,bpx,wfx,0,  l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 
le£=ap(ind)-0.57; 

if  lef<0  ,lef=0;  end 

ap=linspace(lef,  ap(ind)+0 .57,8); 
%  four  more  time  to  find  alpha 

x3=0; 

fori=l:64 

x3=x3+... 

solalph(,xx00,,0,2*pi,2,bpx,wfx,0,l,2,bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+... 

solalph('xxir,0,2*pi,2,bpx,wfx,0,l,2,bpx,wix,ap,b,CT,i,I,k,m,RPTN); 
end 

[val,ind]=min(x3); 
ap=ap(ind); 
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%  after  find  optimal  alpha  use  double  integration 
qq3=0; 

for  i=  1:64 

qq3=qq3+(  dbgquadm('xx00',0,2*pi,2,bpx,  wfx,0, 1,2,... 

bpx,wfx,ap,b,CT,i,I,k,m,RPTN)+. . . 

dbgquadm('xx  1 1  ',0,2*pi,2,bpx,  wfx,0, 1,2,... 

bpx,wfx,ap,b,CT,i,I,k,m,RPTN)  )/(256*pi); 
end 

BER=[BER  qq3];  thresh2=[thresh2  ap]; 

if  qq3<10A(-15)         %for  save  time  only  interesting  in  10A(-15) 
ii=find(ppp=RPTN); 

BER(ii+l:lenl)=5*10A(-116)*ones(  l,length(ii+l:lenl) ); 
thresh2(ii+l  :lenl)=(ap+2)*ones(l,length(ii+l  :lenl) ); 
break 
end 
end 

pe=[pe;BER];  thresh l=[threshl;thresh2]; 
end 

pp=[pp;pe];  thresh=[thresh;threshl]; 
end 

time2=toc; 
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%  worse  case  form  appendix  setting  optimal  threshold  equal  (x0+xl)/2 

M=4; 

I=linspace(0,6, 101);  %  set  I  value  in  linspceinteger  to  find  minimum  I 

RPTN_DB=0:0.2:20;  %  x-axis  dB  range 

RPTN=10A(0. 1  *RPTN_DB);    %  change  to  ratio 

%  single  channel 

BER0=0. 5  *eifc(RPTN/8A0. 5); 

%find  optimal  alpha 

pp=[]; 

forCT=[5  10  15  20] 

x3=0; 

fork=l:0.5*M 

x3=x3+2./(l+(2*pi*k*I/CT)A2); 

end 

x0=(l-exp(-2*CT))/(2*CT)+x3*(l+2*(l-exp(-CT))/CT); 
xl=l-2*(l-exp(-CT))/CT+(l-exp(-2*CT))/(2*CT); 
[val,ind]=min(abs(xO-xl)); 
I=ceil(I(ind)); 
I=(I(ind));  %  our  minimum  I  value 

%  set  different  I  and  to  find  the  BER 
ss=I+l:I+10; 
BER=[]; 
for  I=ss 
x3=0; 
fork=l:0.5*M 

x3=x3+2/(l+(2*pi*k*I/CT)A2); 
end 

x0=(l-exp(-2*CT))/(2*CT)+x3*(l+2*(l-exp(-CT))/CT); 
xl=l-2*(l-exp(-CT))/CT+(l-exp(-2*CT))/(2*CT); 
BER=[BER;0.5*erfc(  RPTN*(xl-x0)/8A0.5  )  ]; 
end 

pp=[pp;BER]; 
end 

semilogy(RPTN_DB,BERO.,--,,RPTN_DB,pp(:,:)) 

axis([10  19  10A(-15)  1]) 
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